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ABSTRACT 

It has long been expected that some massive stars produce stellar mass black holes (BHs) upon death. Un- 
fortunately, the observational signature of such events has been unclear It has even been suggested that the 
result may be an "unnova," in which the formation of a BH is marked by the disappearance of a star rather than 
an electromagnetic outburst. I argue that when the progenitor is a red supergiant, evidence for BH creation 
may instead be a w 3- lOday optical transient with a peak luminosity of sa 10'*°- 10'" ergs"', a temperature of 
~ lO'^K, slow ejection speeds of k, 200 kms"', and a spectrum devoid of the nucleosynthetic products associ- 
ated with explosive burning. This signal is the breakout of a shock generated by the hydrodynamic response 
of a massive stellar envelope when the protoneutron star loses ^ few x 0.1 M0 to neutrino emission prior to 
collapse to a BH. Current and future wide-field, high-cadence optical surveys make this an ideal time to dis- 
cover and study these events. Motivated by the unique parameter space probed by this scenario, I discuss more 
broadly the range of properties expected for shock breakout flashes, with emphasis on progenitors with large 
radii and/or small shock energies. This may have application in a wider diversity of explosive events, from pair 
instability supernovae to newly discovered but yet to be understood transients. 
Subject headings: black hole physics — stars: evolution — stars: transients — supernovae: general 



1. INTRODUCTION 

It is currently unknown what fraction of massive stars pro- 
duce black holes (BHs) rather than neutron stars (NSs), what 
the channels for BH formation are, and what corresponding 
observational signatures are expected. There is strong evi- 
dence for stellar mass BHs from X-ray binaries throughout 
our galaxy (Remillard & Mc Clintock 200 6), so it is clear BHs 
must be a possible endpoint of stellar evolution. Pre-explosive 
imaging of core-collapse su pernovae (SNe) sug gests progen- 
itor masses < 17-20Mq (ISmartt et all 120091) for standai'd 
Type II-P SNe, which are thought to produce NSs. Assum- 
ing a Salpeter initial mass function, this implies that an upper 
limit of ^ 30-35% of massive stars above %Mq fail to lead 
to a successful SNe, and perhaps this number is related to 
the fraction of BHs produced. Such a straightforward com- 
parison is complic ated by the impact of binary interactions 
JSmith et al.ll2011h . which are expected to dom inate the evo- 
lution of most massive stars dSana et al.ll2012l) . A collapsar 
and gamma-ray bur st likely accompanie s some instances of 
BH formation (e.g.. lMacFadven & Woo slev 1999), but these 
are too rare to explain most BHs and are confined to certain 
environments (IStaneket al J 120061: iModjaz et ani2008h . It is 
possible that the signature of BH formation is in fact the dis- 
appearance of a m assive star, or "unnova ," rather than an ac- 
tual SN-like event (iKochanek et alj|2008h . 

On the theoretical side t here is also m u ch uncertaint y . Work 
bv iTimmes et al.l (IT996h. iFWe^ (IT991. iHepr et all (l200l 



lEIdridge & Touli(l2004.IZhang et al.l(l2008h . b'Connor & Otd 
(l2011h . and lUgliano et al.l (120121) attempts to connect the out- 
comes of stellar collapse to the progenitor zero-age main se- 
quence (ZAMS) mass and metallicity. At solar metallicities, 
such models predict BH formation for roughly ~ 10-25% 
of progenitors, but this depends sensitively on many un- 
certain factors such as the treatment of mass loss. Neu- 
trino emission may be a way of inferring co llapse t o a BH 
JBurr ows 1986; Baun igarte et al.l 119961 : lUebendorf er et al.l 
120041: ISumivoshietal.1120071) . but this wifl only be detected 



for especially nearb y events. Rotation may assist in producing 
aSN-Hke signature (I Wooslev & Heged20r2l:lDexter & KasenI 
I2OI2I) and this may enhance gravitational wave emission 
JPiro & Thrane 2012) . but it is not clear in what fraction of 
events rotation is important. If indeed the collapse proceeds 
as a simple ir nplosion, then t he sta r is not expected to brighten 
significantly (IShapiroll 1 989L [l996h . and an unnova-like signa- 
ture is again implied. 

It is important to remember that massive stars which hydro- 
statically for m degene rate iron cores never d irectly collapse 
to BHs ('e.g.. lBurrowsl[T98llOTonnor & OttllMll) . BH for- 
mation is always preceded by a protoneutron star phase with 
abundant emission of neutrinos ('Burrow sll988llBeacom et al.l 
I2OOII) and gravitational waves (Ott 200^ until the protoneu- 
tron star contracts within i ts event horizon. In a somewhat 
forgotten theoretical study, iNadezhiril (119801) focused on the 
impact of the loss of ^ few x O.IM0 over several seconds 
from this neutrino emission. The star reacts as if the grav- 
itational potential of the core has abruptly changed and ex- 
pands in response. This develops into a shock propagating 
into the dying star's envelope. Although the shock's energy of 
E r^ lO''^- 10'*** erg is relatively small in comparison to typical 
core-collapse SNe, it can nevertheless eject the loosely bound 
envelope of a red super giant. This idea was recently investi- 
gated in more detail by iLovegrove & Woosleyl (120131) . Their 
study focused on the SN-like transient that results. Their gen- 
eral conclusion was that such an event would be dominated 
by recombination of hydrogen, similar to Type II-P SNe. This 
would produce a plateau-like light curve lasting a year or more 
with a luminosity of ^ few x 10^''erg s~' and a cool tempera- 
ture of < 4000 K. So although this may be a more promising 
counterpart to BH formation than an unnova, it would still be 
challenging to identify with current observational capabilities. 

Motivated by these previous studies, I investigate in more 
detail the shock breakout expected from this mechanism. This 
demonstrates that the breakout emission may be the most 
promising signature of BH formation from a red supergiant. 
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In ^ I estimate the main properties of the shock breakout, 
and discuss more broadly the emission properties at large radii 
and/or low energies. I conclude in ^with a summary of my 
results and a discussion of future work. 

2. SHOCK BREAKOUT ESTIMATES 

I begin by briefly reviewing the physics of shock breakout. 
For a more detailed background, the interested reader should 
refer to the abundant literature on this topic (including, but not 
limite d to llmshennik & NadezhinllI989l: iMatzner & McKee 
19991 hereafter MM99: iNakar & Sarill2010l I20I21 iKatz et al' 
20I2I) . The main difference in this present work is the rela- 
tively low energy of the shock. 

The envelope is approximated as a polytrope with index 
n given by po = P\{R*/r- 1)", where pi is the half-radius 
density, and I focus on n = 3/2, as is appropriate for a con- 
vective envelope. As is common practice, I scale p\ by 
the characteristic density of the ejecta /?» = M^JR\, where 
Mej is the mass of the ejecta, since the ratio pi/p* is typ- 
ically of order unitjQ- For the very outer parts of the star, 
I use the variable x = I- r/R^, and set po ~ P\x", where 
X <^ 1. The shock propagates through the envelope with 
speed y, = T(E/my^-(m/poPf (MM99), where E is the 
shock energy, r is the radial coordinate, m is the mass in- 
terior to r, and po is the density profile of the envelope 
prior to expansion due t o the shock. Using the s elf-similar, 
pla nar soluti ons of Gande l'Man & Frank-Kamene tskii ( 1956) 
and lSakurail(IT960l) . I adopt the values of /3 = 0.19 and 0.23 
when the shock is radiation-pressure and gas-pressure dom- 
inated, respectively. The constant F can be estimated from 
self-similar blastwave solutions (MM99), for which I use 
F = 0.794. 

Breakout happens when a shock gets too close to the stellar 
surface, and photons diffuse out into space at an optical depth 
r « c/\'s. Setting r = npxR^x"^^ /{n+ 1), breakout occurs at 
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where v* = {E/Mq)^l-. ScaHng this to typical values for a low 
energy shock in a massive progenitor. 
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where ko34 = K./03Acm^g-\ £48 = £'/10'*'^erg, Mio = 
Mej/1OM0, and/?iooo=-R*/1000/?o. 

The amount of energy in the radiation field associated with 
the shock depends on whether or not the shock is domi- 
nated by radiation pressure. For the radiation-dominated 
case, the shock jump condition is aT^ jT) = 6pQV^/7, where 
a is the radiation constant. Setting the energy to be Ei-^d « 
ATrRlxboiaT"^ /3), where Xbo is the depth at which r = Tbo, 
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If the shock is gas-dominated then pohgT / finip = 3poVj/16, 



' See the Appendix of lCalzavara & Matzneil2004l for derivations of values 
for pi/ Pa, for different stellar structures. 



where ^b is Boltzmann's constant. In this cas^ 
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The actually shock breakout energy is roughly given by ^bo ~ 
min(£'rad,£'gas). A key point is the strong scalings of Egas oc 
£■4-4^ which suppresses shock breakout for small E. For this 
reason, all equations in this paper with numerical factors as- 
sume the radiation-dominated regime unless otherwise noted. 
The observed luminosity is determined by the timescale 
over which this energy i s emitted. The two domin ant effects 
are the light-travel time (lEnsman & Burrowslll992b . 



fit «7;*/c = 39 7^1000 min, 
and the diffusion time (MM99; Piro & Nakar 2012), 
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The timescale for the breakout emission is fbo ~ rnax(fit,fdiff)- 
Note that fdiff > h\ when 
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so that for the low energy shock investigated here fbo ~ fdiff- 
If fdiff is too long, then the breakout is suppressed by cooling 
from adiabatic expansion. This occurs on a timescale. 
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which includes a factor of 2 from geometric effects and pres- 
sure gradients (MM99). 

An additional issue is thermalization (see the discussion 
in Nakar & Sari 2010) . When thermal emission is achieved 
at temperature T , the number density of photons is «ph ~ 
aT^ /3kB- The rate of photon production from free-free emis- 
sion is «ph « 3.5 x 10"'*p^r"'''^s"' cm"-'. Therefore thermal- 
ization is expected for times later than roughly 
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Since ftherm ^ fit, fdiff at low energies, the shock breakout is in 
the thermal regime. This is in contrast to typical SNe with an 
energy of -- 10^' erg where fbo ~ fit and ftherm > fbo- 

In Figure [1] I highlight the main physical processes that 
determine the shock breakout emission as a function of /?* 
and E. The upper panel is for an ejecta mass of Mgj = 3Mq 
and the lower panel is for Mej = IOM0. The diagonal solid 
line divides where the breakout timescale is determined by f^ 
or fdiff. The dotted line shows where the emission is expected 
to transition to being non-thermal (blue, lightly shaded). The 

^ Note that although I denote this energy with a subscript "gas," this energy 
is still associated with the radiation field. This merely identifies that gas 
pressure is dominant in this regime. 
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Fig. 1. — Main physical conditions that deteiTnine shock breakout as a 
function of R„ and E. Top and bottom panels are ejecta masses of 3Mq 
and IOMq, respectively. The solid line divides regions where the duration 
fbo is determined by either the light-crossing time ?i, or the thermal diffu- 
sion time ^(jiff. Above the dotted line is the non-thermal regime (blue, lightly 
shaded). Below the dashed hne the shock becomes gas-pressure dominated 
(red, darkly shaded). In the far right bottom comer, fexp < ?bo and the breakout 
is suppressed by adiabatic expansion. 

dashed Hne shows where the shock starts to be gas pressure 
dominated (red, darkly shaded). 

In Figure |2l I show how these different physical conditions 
translate into observed luminosities and durations. Solid lines 
denote constant breakout luminosity, which is set as 

(11) 



Lbo = min(£'i-ad,£'gas)/max(fit,fdiff), 



where 1 add an interpolation in each of the min and max 
functions to smooth the transitions. Furthermore, 1 include 
a fac tor of \\ + {tho I W tiS^Y^ to account for adiabatic expan- 
sion (iPiro et al.ll2010l) . where 7 is the adiabatic exponent. The 
sohd lines are labeled by the x, where Lbo = lO'ergs"'. The 
breakout luminosity decreases rapidly once gas pressure dom- 
inates, and 1 shade regions where Lbo < 10^^ ergs"'. Dotted 
lines donate constant breakout durations fbo = inax(rit,fdiff)- 
Red crosses compare co re-collapse of a blue supergiant 
jEnsman & Burrowslll992h with BH formation from a red su- 
pergiant, demonstrating the different regimes these two cases 
occupy. 

Since it may be generally useful, a wider range of R^ and 
E are plotted in Figures [T] and |2]than just what applies to the 
BH formation case that is the focus here. For example, for a 
pair instabihty SN with R^ - IQQORq and £ - 10^^- 10^^ erg, 
the expectation from Figure |2]is that Lbo '^ IC'^ergs"' with 
fbo ~ hrs . This is roughly co nsistent with more detailed calcu- 
lationfl jKasen et al.ll201ll) . Similarly, if the luminosity and 
duration of a purported breakout flash is observed, these re- 
sults can be used to infer R^, and E, potentially even providing 
constraints on transient events that are difficult to classify. 

^ The progenitors of pair instability SNe are considerably more massive, 
but this is easily corrected for by using the scalings provided in this work. 
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Fig. 2. — Same as Figure [T] this time plotting contours of constant break- 
out luminosity Lbo (solid lines) and constant duration /bo (green, dotted lines). 
The contours are labeled with their values, with the contours of constant lu- 
minosity labeled by x where Lbo = 10' ergs"'. Purple, shaded regions de- 
note where Lbo < lO^'ergs"'. The decrease in Lbo at especially large R, 
and E (the upper-left and lower-right comers) is due to adiabatic expansion. 
The up per red cross shows a characteristic model from lEnsman & Burrowsl 
09921) . which is appropiiate for core-collapse of a blue supergiant. The lower 
red cross con'esponds to BH formation of a 15 Mq ZAMS red supergiant 
jLovegrove & Woosle\l20Ol maximum mass loss model). 



Especially relevant to the situation of BH formation is the 
regime where Lbo ~ E^i and fbo ~ fdiff , which results in 
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This applies to the entire lower-right triangle of parameter 
space in each of the panels of Figure [T] (although note that 
adiabatic expansion can cause the luminosity to be some- 
what lower than this). The observed shock breakout temper- 
ature is in the thermal regime and thus can be estimated as 

Tobs ~ T^jo/t^jo , where Lbo is the temperature of the plasma at 
the depth of the breakout, resulting in 
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Since 7bo > Lobs, an electron scattering opacity is sufficient to 
estimate the diffusion properties. At the surface, where the 
temperature is close to Lobs, other opacity effects may impact 
the observed spectra, which should be modeled in the future. 
Subsequent to the breakout emission, there is a plateau- 
like phase due to hydr o gen re combination as investigated by 
iLovegrove & Woosleyl (120 13h . Using the analytic fits to the 
numerical models of Type II-P SNe by iKasen & Woosleyl 
(l2009h . I estimate that this phase has a luminosity 
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with a duration of 



fplal 



i 420 '",!°"'' days, 

^48 



(15) 



which roughly matches iLove grove & Wooslevl(l2013l) . Com- 
parison between Lbo and Lpiat shows that breakout emission 
will be more conducive to detection by surveys. 

An important constraint on observing the breakout from BH 
formation is that the shock must have sufficient energy to pro- 
duce an observable signal. Estimating the breakout duration 
in the gas-dominated regime, 
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Setting this as > lO^'^ergs ', the shock energy must obey 
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iLovegrove & WooslevI (120 13h find values around this range, 
both above and below, and with generally more energy for 
a 15 Mq ZAMS star in comparison to 25 Mq. This suggests 
that the breakout may not be detectable in all cases. An im- 
portant focus for future simulations of collapsing stars is to 
better explore the full range of shock energies possible. 

3. CONCLUSIONS AND DISCUSSION 

I have investigated shock breakout during the collapse of 
a red supergiant producing a BH. The shock is generated 
by the hydrodynamic response of ^ few x O.IMq being car- 
ried away by neutrinos in the time before BH formation. 
This breakout flash has many distinctive characteristics that 
will help distin guish it from other po te ntial short timescale 
transients (e.g.. iMetzger et al.ll2009allbl: iDarbhaet all 1201 Ot 
iMetzger et al.ll2010tlShen et al.ll2010l) . which I summarize as 
follows. 

1. The breakout flash has a luminosity Lbo ~ 10^"- 
lO'^ergs"', as given by equation (fT2l i. 

2. The flash duration is fbo ~ 3- lOdays with fbo oc /?J '^ 
because it is set by fdiif (rather than fic as in many SNe). 

3. The breakout is in the thermal regime with an observed 
temperature of Tobs ~ 10,000K. 



4. The velocities should be relatively low with v,nax ~ 
200kms~', as given by equation (|9]l. 

5. The star likely retained most of its hydrogen envelope 
to have a sufficiently large radius for a bright and long 
lasting breakout, and this would be seen in the spectra. 



6. The spectrum should be devoid of nucleosynthetic 
products from explosive burning. 

The source would peak in the ultraviolet with an absolute 
magnitude of roughly -14.5. The spectrum will be bright in 
blue and visual wavebands, making it well-suited for detec- 
tion by wide-field, transient surveys like the Palorn ar Tran- 
sient Factory (PTF:l Rau et"aI1l2009l: ITaw et al.ll2009h and the 
Panoramic Sur vey Telescope and Rapid Response System 
(Pan-STARRS: iKaiser et"al]l2002h . Future theoretical work 
should better quantify the time-dependent temperature and lu- 
minosity during this phase. If the progenitor is a blue super- 
giant or a Wolf-Rayet star rather than a red supergiant, then 
the smaller radius causes the shock breakout duration to be 
significantly shorter than what I summarize here. Neverthe- 
less, the main properties can be estimated using Figure |2] 

A critical issue I have explored, which typically does not 
occur for normal shock breakouts from SNe, is the impact 
of gas pressure and adiabatic expansion at low shock ener- 
gies. It is shown that if the shock is sufficiently low en- 
ergy, as quantified by equation ( fTSl l. then the breakout will 
not be observable. A critical question for future theoretical 
studies will therefore be to explore what exactly is the ex- 
pected range of energies for this shock. This may depend 
on many factors, including the progenitor mass, neutron star 
equation of state, treatment of neutrino emission, and dura- 
tion of neutrino emission from the protoneutron star p rior to 
collapse to a BH. For example. iLovegrove & Woosley (' 2013h 
find that a NS equation of state that favors a massive maxi- 
mum mass (« 2.5Mq) results in increased neutrino emission 
and a stronger shock. This shows how the detection and study 
of the signal described here could assist in addressing other 
fundamental questions in physics and astrophysics. 
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